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The role of surface collective excitations in the electron relaxation in small 
metal particles is studied. It is shown that the dynamically screened electron- 
electron interaction in a nanoparticle contains a size-dependent correction in- 
duced by the surface. This leads to new channels of quasiparticle scattering 
accompanied by the emission of surface collective excitations. In noble-metal 
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particles, the dipole collective excitations (surface plasmons) mediate a res- 
onant scattering of (i-holes to the conduction band. The role of this effect 



C — \ in the ultrafast optical dynamics of small nanoparticles is studied. With 



decreasing nanoparticle size, it leads to a drastic change in the differential 
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^> ■ absorption lineshape and a strong frequency dependence of the relaxation 

near the surface plasmon resonance. The experimental implications of these 

-j*^ . results in ultrafast pump-probe spectroscopy are addressed. We also discuss 

the size-dependence of conduction electron scattering rates. 
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Surface col ective excitations play an important role in the absorption of light by metal 
nanoparticles.B~i In large particles with sizes comparable to the wave-length of light A (but 
smaller than the bulk mean free path), the lineshape of the surface plasmon (SP) resonance 
determined by the electromagnetic effects.0 On the other hand, in small nanoparticles 
with radii R <C A, the absorption spectrum is governed by quantum confinement effects. For 
example, the momentum non-conservation due to the confining potential leads to the Landau 
damping of the SP and to a resonance linewidth inversely proportional to the nanoparticle 
sizefli Confinement also changes non-linear optical properties of nanoparticles: a size- 
dependent enhancement of the third-order susceptibilities, caused by the elastic surface 
scattering of single-particle excitations, has been reportedJlrEI 

Recently, extensive experimental studies of the electron relaxation in nanoparticles have 
been performed using ultrafast pump-probe spectroscopy.0^0 Unlike in semiconductors, 
the dephasing processes in metals are very fast, and nonequilibrium populations of optically 
excited electrons and holes are formed within several femtoseconds. These thermalize into 
the hot Fermi-Dirac distribution within several hundreds of femtoseconds, mainly due to 
e-e and h-h scattering]!!^ Since the electron heat capacity is much smaller than that of the 
lattice, a high electron temperature can be reached during less than 1 ps time scales, i.e., be- 
fore any significant energy transfer to the phonon bath occurs. Note that in nanometer-sized 
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metal particles and nanoshels, the electron-phonon coupling is weaker than in the bulk.E30 
During this stage, the SP resonance has been observed to undergo a time-dependent spec- 
tral broadeningB0 Subsequently, the electron and phonon baths equilibrate through the 
electron-phonon interactions over time intervals of a few picoseconds. During this incoherent 
stage, the hot electron distribution can be characterized by a time-dependent temperature. 

The many-body correlations play an important role in the transient changes of the 
absorption spectrum.BHi] For example, in undoped materials, four-particle correlations and 
exciton-exciton interactions were shown to play a dominant role in the optical response 
for specific sequences of the optical pulses.0c! Correlation effects also play an important 
role in doped systemsilE In nanoparticles, it has been shown that in order to explain the 
differential absorption lineshape, it is essential to take into account the energy-dependent e- 
e scattering of the optically-excited carriers near the Fermi surfaced Furthermore, despite 
the similarities to the bulk-like behavior, observed, e.g., in metal films, certain aspects of the 
optical dynamics in nanoparticles are significantly different For example, experimental 
studies of small Cu nanoparticles revealed that the relaxation times of the pump-probe signal 
depend strongly on the probe frequency: the relaxation was considerably slower at the SP 
resonance. @L3 This and other observations suggest that collective surface excitations play 
an important role in the electron dynamics in small metal particles. 

In this paper we address the role of surface collective excitations in the electron re- 
laxation in small metal particles. We show that the dynamically screened e-e interaction 
contains a correction originating from the surface collective modes excited by an electron 
in nanoparticle. This opens up new quasiparticle scattering channels mediated by surface 
collective modes. We derive the corresponding scattering rates, which depend strongly on 
the nanoparticle size. The scattering rate of a conduction electron increases with energy, in 
contrast to the bulk-plasmon mediated scattering. In noble metal particles, we study the 
SP-mediated scattering of a <i-hole into the conduction band. The scattering rate of this 
process depends strongly on temperature, and exhibits a peak as a function of energy due 
to the restricted phase space available for interband scattering. We show that this effect 
manifests itself in the ultrafast nonlinear optical dynamics of nanometer-sized particles. We 
perform self-consistent calculations of the temporal evolution of absorption spectrum in the 
presence of the pump pulse. For large sizes, the absorption peak exhibits a red shift at short 
time delays. We show that with decreasing size, the SP resonance developes a blue shift 
due to the SP-mediated interband scattering of the <i-hole. We also find that the relaxation 
times of the pump-probe signal depend strongly on the probe frequency, in agreement with 
recent experiments. 

The paper is organized as follows. In Section |TJ we review the relevant basic relations 
of the linear response theory for nanoparticles. In Section |l| the dynamically screened 
Coulomb potential in a nanoparticle is derived. In Section |IV] we calculate the SP-mediated 
interband scattering rate of a d-band hole. In Section [V] we incorporate this effect in the 
calculation of the absorption coefficient. In Section |VI| we present our numerical results and 
discuss their experimental implications on the size and frequency dependence of the time- 
resolved pump-probe signal. In Section [VII| we study the quasiparticle scattering in the 
conduction band mediated by surface collective modes. Section |V111| concludes the paper. 
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II. BASIC RELATIONS 



In this section we recall the basic relations regarding the linear absorption by metal 
nanoparticles embedded in a medium with dielectric constant e m . We will focus primarily 
on noble metal particles containing several hundreds of atoms; in this case, the confinement 
affects the extended electronic states after the bulk lattice structure has been established. 
If the particles radii are small, C A, only dipole surface modes can be optically excited 
and non-local effects can be neglected. In this case the optical properties of this system are 
determined by the dielectric functiori0 

e co i(w) = e m + 3pe m — — — — , (1) 
€{w) + 2e m 

where e(a>) = e'(uj) + ie"(u) is the dielectric function of a metal particle and p <C 1 is 
the volume fraction occupied by nanoparticles in the colloid. Since the (i-electrons play an 
important role in the optical properties of noble metals, the dielectric function e(u) includes 
the interband contribution 6d(u>). For p <C 1, the absorption coefficient of such a system is 
proportional to that of a single particle and is given bytl 

«M = -9pe^ 2 -Im^-, (2) 

C € S {UJ) 

where 

e e (u) = e d (u) -u 2 /uj(uj + ij s ) + 2e m , (3) 

plays the role of an effective dielectric function of a particle in the medium. Its zero, 
e' s (uj s ) = 0, determines the frequency of the SP, u> s . In Eq. (||), u p is the bulk plasmon 
frequency of the conduction electrons, and the width 7 S characterizes the SP damping. The 
semiclassical result Eqs. (||) and (0) applies to nanoparticles with radii R 3> q~ F , where q TF 
is the Thomas-Fermi screening wave-vector ~ 1 A in noble metals). In this case, the 
electron density deviates from its classical shape only within a surface layer occupying a small 
fraction of the total volume.^ Quantum mechanical corrections, arising from the discrete 
energy spectrum, lead to a width 7 S ~ v F /R, where v F = k F /m is the Fermi velocity.^! 
Even though 7 s /co> s ~ (g TF i?) _1 <C 1, this damping mechanism dominates over others, e.g., 
due to phonons, for sizes R < 10 nm. In small clusters, containing several dozens of atoms, 
the semiclassical approximation breaks down and density functional or ab initio methods 
should be used.S^i 

It should be noted that, in contrast to surface collective excitations, the e-e scattering 
is not sensitive to the nanoparticle size as long as the condition q TF R S> 1 holds.E^ Indeed, 
for such sizes, the static screening is essentially bulk-like. At the same time, the energy 
dependence of the bulk e-e scattering rate,0 7 e oc (E—E F ) 2 , with E F being the Fermi energy, 
comes from the phase-space restriction due to the momentum conservation, and involves the 
exchange of typical momenta q ~ q TF . If the size-induced momentum uncertainty 5q ~ 
is much smaller thano TF , the e-e scattering rate in a nanoparticle is not significantly affected 
by the confinement.^ 
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III. PLASMON POLE APPROXIMATION IN SMALL METAL PARTICLES 



In this section, we study the effect of the surface collective excitations on the e-e inter- 
actions in a spherical metal particle. To find the dynamically screened Coulomb potential, 
we generalize the method previously developed for calculations of local field corrections to 
the optical fields.S The potential U(lu; r, r') at point r arising from an electron at point r' 
is determined by the equation!! 

U (uj; r, r') = u(r — r') + / dr 1 dr 2 u(r - n)II(a;; ri, r 2 )U(u; r 2 , r'), (4) 

where u(r — r') = e 2 |r — r'| _1 is the unscreened Coulomb potential and n(c<j; i - !, r 2 ) is the 
polarization operator. There are three contributions to II, arising from the polarization of 
the conduction electrons, the (/-electrons, and the medium surrounding the nanoparticles: 
II = Il c + ILj + II m . It is useful to rewrite Eq. (H) in the "classical" form 

V- (E + 4?rP) = 47re 2 d(r-r'), (5) 

where E(c<j; r, r') = —VU{oj; r, r') is the screened Coulomb field and P = P c + + P m is 
the electric polarization vector, related to the potential U as 

VP{u; r, r') = -e 2 J dr^u; r, n)l/(a;; r u r'). (6) 

In the random phase approximation, the intraband polarization operator is given by 

n c (a;;r,rO = £ ^ ( ^"/^^ ^(r)V»g(r)V»r(r , )^(r , ) > (7) 

where E c a and ip^, are the single-electron eigenenergies and eigenfunctions in the nanopar- 
ticle, and f(E) is the Fermi-Dirac distribution (we set h = 1). Since we are interested in 
frequencies much larger than the single-particle level spacing, H c (uj) can be expanded in 
terms of 1/uj. For the real part, W c (uj), we obtain in the leading ordered 

U' c (u; r, n) = l —V[n c (r)V6(r - n)], (8) 

muj z 

where n c (r) is the conduction electron density. In the following we assume, for simplicity, 
a step density profile, n c (r) = n c 6(R — r), where n c is the average density. The leading 
contribution to the imaginary part, H"(u), is proportional to u~ 3 , so that U."(cj) <C U' c (uj). 
By using Eqs. (|S|) and (^), one obtains a familiar expression for P c at high frequencies, 

P c (u;- r, r') = f^Vtf (w; r, r') = 6(R - r)x c (w)E(w; r, r'), (9) 

where Xc{w) = —e 2 n c /muj 2 is the conduction electron susceptibility. Note that, for a step 
density profile, P c vanishes outside the particle. The <i-band and dielectric medium contri- 
butions to P are also given by similar relations, 

P d (oo; r, r') = G(R - r) Xd (oj)E(uj; r, r'), (10) 
P m (u- r, r') = 9(r - i?)x m E(^; r, r'), (11) 
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where Xi — ( e t — l)/4vr, i = d,m are the corresponding susceptibilities and the step functions 
account for the boundary conditions.il Using Eqs. (Pl)- (|Tl"l) , one can write a closed equation 
for U(uj; r, r'). Using Eq. (|6]), the second term of Eq. (H) can be presented as — e~ 2 J dviu(r — 
r x ) V • P(uj; ri, r'). Substituting the above expressions for P, we then obtain after integrating 
by parts 

e(uj)U(uj; r, r') = / dnVx— r • Vi [9(R - r) X (uj) + 9(r - R) Xm ] U(lu; n, r') 

|r — r I J |r — ri| 

+i [ G?r l0 ?r 2 — ^ i U'^cu;r 1 ,r 2 )U(oj;r 2 ,r'), (12) 

J \r — Ti\ 

with 

e(w) = 1 + 4<k X {u) = e d (u) - uj 2 p /uj 2 , (13) 

uj 2 = 47re 2 n c /m being the plasmon frequency in the conduction band. The last term in 
the rhs of Eq. (0), proportional to IT" (a>), can be regarded as a small correction. To solve 
Eq. fljjD , we first eliminate the angular dependence by expanding U(u;r, r') in spherical 
harmonics, Y LM (r), with coefficients UiM(^', r , r ')- Using the corresponding expansion of 

|r — r'| _1 with coefficients QLAi( r , r ') — 2L+i r L lr ' L (^ or r > r ')> we ^ ne following 
equation for ULM(u;r,r'): 

L -\- 1 f t \^ 

e(u)U LM (uj; r, r') = Q LM (r, r') + 4?r [x(u;) - Xm] 2 L + I \RJ Ulm ^' R > r "> 

+ ie> £ / ( ir 1 * 2 r?r5Q IM (r,r 1 )n' iM , 1 , M ,( W ;r 1 , r2 )^ M ,( W ;r 2 , r 0, (14) 

L'M' 7 

where 

n LAf,L'M'( w ; r i ; r 2) = / rfrirfr2y2M(fi)n"(^;ri,r 2 )yL/Af'(r2), (15) 

are the coefficients of the multipole expansion of II" (cu; ri, r 2 ). For U" = 0, the solution of 
Eq. (14) can be presented in the form 

U LM {u>;r,r') = a(uj)e 2 Q LM (r, r') + b(u) ^ - - R2L+l , (16) 

with frequency-dependent coefficients a and b. Since IT"(co>) <C n' c (u;) for relevant frequen- 
cies, the solution of Eq. fll~4| ) in the presence of the last term can be written in the same 
form as Eq. (|I"6D, but with modified a{oj) and b{u). Substituting Eq. ( |i~6"l) into Eq. (fill) , we 
obtain after lengthy algebra in the lowest order in II" 

aH = e" 1 H, b(u) = ei\u)-e-\u), (17) 

where 

= 2L + i e ^ + 2TTi em + ie '^ u) ' (18) 
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is the effective dielectric function, whose zero, e' L (u>L) = 0, determines the frequency of the 
collective surface excitation with angular momentum L,0 

Ul = Le' d (cu L ) + (L + l)e m (19) 

In Eq. ([18|) , e" L (ou) characterizes the damping of the L-pole collective mode by single-particle 
excitations, and is given by 

= {2L - 1)jR2L+1 £ \M™\ 2 [f(K) - f(E c a ,)}5(E c a - E c a , + uj), (20) 

where Mj$ are the matrix elements of r L >£A,/(r). Due to the momentum nonconservation 
in a nanoparticle, the matrix elements are finite, which leads to the size-dependent width 
of the L-pole mode: 

2L + 1 u 3 ., . . 
1l = — (21) 

For uj ~ u>l, one can show that the width, ji ~ v F /R, is independent of uj. Note that, in 
noble metal particles, there is an additional d-electron contribution to the imaginary part 
of ei(u) at frequencies above the onset A of the interband transitions. 

Putting everything together, we arrive at the following expression for the dynamically- 
screened interaction potential in a nanoparticle: 

u ^ - Uj ^f + 5§ ^TTiZR {if™™. (22) 

with e~l l (uj) = e^ 1 (u;) — e~ l (uj). Equation (p2[) , which is the main result of this section, 
represents a generalization of the plasmon pole approximation to spherical particles. The 
two terms in the rhs describe two distinct contributions. The first comes from the usual 
bulk-like screening of the Coulomb potential. The second contribution describes a new 
effective e-e interaction induced by the surface: the potential of an electron inside the 
nanoparticle excites high-frequency surface collective modes, which in turn act as image 
charges that interact with the second electron. It should be emphasized that, unlike in the 
case of the optical fields, the surface-induced dynamical screening of the Coulomb potential 
is size-dependent. 

Note that the excitation energies of the surface collective modes are lower than the bulk 
plasmon energy, also given by Eq. ( |19|) but with e m = 0. This opens up new channels of 
quasiparticle scattering, considered in the next section. 



IV. INTERBAND SCATTERING MEDIATED BY SURFACE PLASMONS 

We now turn to the interband processes in noble metal particles and consider the scat- 
tering of a d-ho\e into the conduction band. We restrict ourselves to the scattering via the 
dipole channel, mediated by the SP. The corresponding surface-induced potential, given by 
the L = 1 term in Eq. (|22"D , has the form 
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^p 2 r • r' 1 

u -^ = -r-*1M- (23) 

With this potential, the ci-hole Matsubara self-energy is given by0 

where d aQ / = (c, a\r\d, a') = (c, a\p\d, a') / im(E^ — E^,) is the interband transition matrix 
element. Since the final state energies in the conduction band are high (in the case of 
interest here, they are close to the Fermi level), the matrix element can be approximated by 
the bulk-like expression (c, a\p\d, a') = 5 aa i(c\p\d) = 5 aa >fi, the corrections due to surface 
scattering being suppressed by a factor of (/c F i?) _1 <C 1. After performing the frequency 
summation, we obtain for ImS^ 

ImE'M- 9eV N(E° a - u) + f(ET a ) 

ImSoM " "^Wp Ira ' (25) 

with E^ = E c a — E^. We see that the scattering rate of a ci-hole with energy E%, r yf l (E^,) = 
ImS^E^), has a strong R~ 3 dependence on the nanoparticle size. 

The important feature the interband SP-mediated scattering rate is its energy depen- 
dence. Since the surface-induced potential, Eq. (p3|), allows for only vertical (dipole) inter- 
band single-particle excitations, the phase space for the scattering of a ci-hole with energy 
E% is restricted to a single final state in the conduction band with energy E c a . As a result, 
the ci-hole scattering rate, ■jf l (E^,), exhibits a peak as the difference between the energies of 
final and initial states, E^ = E c a — E%, approaches the SP frequency u s [see Eq. (|25])]. 

As we show in the next section, the fact that the scattering rate of a ci-hole is dominated 
by the SP resonance, affects strongly the nonlinear optical dynamics in small nanoparticles. 
This is the case, in particular, when the SP frequency, uj*. is close to the onset of interband 
transitions, A, as, e.g., in Cu and Au nanoparticles.liii3Mll3 Consider an e-h pair with 
excitation energy uj close to A. As we discussed, the ci-hole can scatter into the conduction 
band by emitting a SP. According to Eq. (J2^|), for uj ~ lu s , this process will be resonantly 
enhanced. At the same time, the electron can scatter in the conduction band via the usual 
two-quasiparticle process. For uj ~ A, the electron energy is close to Ep, and its scattering 
rate is estimated asll 7 e ~ 1CT 2 eV. Using the bulk value of fi, 2/j 2 /m ~ 1 eV near the 
L-point,0 we find that 7^ exceeds 7 e for R < 2.5 nm. In fact, one would expect that, in 
nanoparticles, fi is larger than in the bulk due to the localization of the conduction electron 
wave-functions.0 



V. SURFACE PLASMON OPTICAL DYNAMICS 

In this section, we study the effect of the SP-mediated interband scattering on the non- 
linear optical dynamics in noble metal nanoparticles. When the hot electron distribution 
has already thermalized and the electron gas is cooling to the lattice, the transient response 
of a nanoparticle can be described by the time-dependent absorption coefficient a(u,t), 
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given by Eq. (|2]) with time-dependent temperature. In noble-metal particles, the temper- 
ature dependence of a originates from two different sources. First is the phonon-induced 
correction to j a , which is proportional to the lattice temperature 7](t). As mentioned in the 
Introduction, for small nanoparticles this effect is relatively weak. Second, near the onset of 
the interband transitions, A, the absorption coefficient depends on the electron temperature 
T(t) via the interband dielectric function €d(oo) [see Eqs. (§) and @]. In fact, in Cu or Au 
nanoparticles, uj s can be tuned close to A, so the SP damping by interband e-h excitations 
leads to an additional broadening of the absorption peak.0 In this case, the temperature 
dependence of e d {uj) dominates the pump-probe dynamics. Below we show that, near the 
SP resonance, both the temperature and frequency dependence of e^w) = 1 + Anxdiu) are 
strongly affected by the SP-mediated interband scattering. 

For non-interactingelectrons, the interband susceptibility, Xdii^>) = Xd(iu) + Xd(-i^), 



has the standard formKH 

M**) = -E^^p^E^^O^K + io;), (26) 

where G^{iuj') is the Green function of a ci-electron. Since the ci-band is fully occupied, the 
only allowed SP-mediated interband scattering is that of the <i-hole. We assume here, for 
simplicity, a dispersionless ei-band with energy E d . Substituting G d (ioo') = [iu' — E d + Ep — 
T^iioj')}^ 1 , with £^(iu;) given by Eq. (fHl) , and performing the frequency summation, we 
obtain 

_ eV r dE<g(E c ) f(E c )~l 

Xd{ ' m 2 J (E cd ) 2 u- E^ + iYh^^E^Y 1 ' 

where g(E c ) is the density of states of conduction electrons. Here 7^(^, E c ) = lmH d (E c — to) 
is the scattering rate of a ci-hole with energy E c — u, for which we obtain from Eq. (|25p, 



where we neglected N(uj) for frequencies uj ~ u s ^> ksT. Remarkably, 7^(0;, E c ) exhibits 
a sharp peak as a function of the frequency of the probe optical field. The reason for this 
is that the scattering rate of a d-hole with energy E depends explicitly on the difference 
between the final and initial states, E c — E, as discussed in the previous section: therefore, 
for a ci-hole with energy E = E c — 00, the dependence on the final state energy, E c , cancels 
out in e s (E c — E) [see Eq. (pop]. This implies that the optically-excited <i-hole experiences 
a resonant scattering into the conduction band as the probe frequency u approaches the SP 
frequency. It is important to note that ■jf l (ui,E c ) is, in fact, proportional to the absorption 
coefficient a(u>) [see Eq. (§)]. Therefore, the calculation of the absorption spectrum is a 
self-consistent problem defined by Eqs. (fj), (|3|), (|27|), and (|28|). 

It should be emphasized that the effect of 7^ on e^(uj) increases with temperature. Indeed, 
the Fermi function in the rhs of Eq. (^) implies that 7^ is small unless E c — Ep < ksT. 
Since the main contribution to x'd(u}) comes from energies E c — Ep ~ cu — A, the ci-hole 
scattering becomes efficient for electron temperatures fc#T > u s — A. As a result, near the 
SP resonance, the time evolution of the differential absorption, governed by the temperature 
dependence of a, becomes strongly size-dependent, as we show in the next section. 
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VI. NUMERICAL RESULTS 



In the numerical calculations below, we adopt the parameters of the experiment of Ref. 
|T2"| , which was performed on R ~ 2.5 nm Cu nanoparticles with SP frequency, uj s ~ 2.22 eV, 
slightly above the onset of the interband transitions, A ~ 2.18 eV. In order to describe the 
time-evolution of the differential absorption spectra, we first need to determine the time- 
dependence of the electron temperature, T(t), due to the relaxation of the electron gas to 
the lattice. For this, we employ a simple two-temperature model, defined by heat equations 
for T(t) and the lattice temperature T}(t): 

C(T)^ = -G(T-T l ), 

C l -^ = G(T-T l ), (29) 

where C{T) = TT and C% are the electron and lattice heat capacities, respectively, and G 
is the electron-phonon coupling.^ The parameter values used here were G = 3.5 x 10 16 
Wm _3 K _1 , T = 70 Jm _3 K -2 , and C\ = 3.5 Jm _3 K _1 . The values of 7 S and /i were extracted 
from the fit to the linear absorption spectrum, and the initial condition for Eq. was 
taken as To = 800 K, the estimated pump-induced hot electron temperature. El We then self- 
consistently calculated the time-dependent absorption coefficient a(u,t), which describes the 
spectrum in the presence of the pump field. The differential transmission is then proportional 
to a r (oj) — a(u,t), where at r (u) is the absorption coefficient at the room temperature. 

In Fig. 1 we show the calculated absorption spectra for different nanoparticle sizes. Fig. 
1(a) shows the spectra at several time delays for R = 5.0 nm; for this size, the SP-mediated 
(i-hole scattering has no effect. With decreasing nanoparticle size, the linear absorption 
spectra are not significantly altered, as can be seen in Figs. 1(b) and (c)]. However, the 
change becomes pronounced at short time delays corresponding to higher temperatures [see 
Figs. 1(b) and (c)]. This effect is clearly seen in the differential transmission spectra, shown 
in Fig. 2, which undergoes qualitative transformation with decreasing size. 

Note that it is necessary to include the intraband e-e scattering, in order to reproduce 
the differential transmission lineshape observed in the experiment. t3 For optically excited 
electron energy close to Ep, this can be achieved by adding the e-e scattering rate0 7 e (-E c ) oc 
[1 - f(E c )][(E c - E F f + (nk B T) 2 ] to Yh in Eq. (g7|). The difference in 7 e (£ c ) for E c below 
and above Ep leads to a lineshape similar to that expected from the combination of red-shift 
and broadening [see Fig. 2(a)]. 

In Figs. 2(b) and (c) we show the differential transmission spectra with decreasing 
nanoparticle size. For R = 2.5 nm, the apparent red-shift is reduced [see Fig. 2(b)]. This 
change can be explained as follows. Since here uo s ~ A, the SP is damped by the interband 
excitations. This broadens the spectra for uo > u s , so that the absorption peak is asymmet- 
ric. The (i-hole scattering with the SP enhances the damping; since the ^-dependence of 7^ 
follows that of a, this effect is larger above the resonance. On the other hand, the efficiency 
of the scattering increases with temperature, as discussed above. Therefore, for short time 
delays, the relative increase in the absorption is larger for u > uj s . With decreasing size, 
the strength of this effect increases further, leading to an apparent blue-shift [see Fig. 2(c)]. 
Such a strong change in the absorption dynamics originates from the R~ 3 dependence of the 



9 



d-ho\e scattering rate; reducing the size by the factor of two results in an enhancement of 
7^ by an order of magnitude. 

In Fig. 3 we show the time evolution of the differential transmission at several frequencies 
close to uj s . It can be seen that the relaxation is slowest at the SP resonance; this char- 
acterizes the robustness of the collective mode, which determines the peak position, versus 
the single-particle excitations, which determine the resonance width. For larger sizes, at 
which 7^ is small, the change in the differential transmission decay rate with frequency is 
smoother above the resonance [see Fig. 3(a)]. This stems from the asymmetric lineshape 
of the absorption peak, mentioned above: the absorption is larger for ui > u s , so that its 
relative change with temperature is weaker. For smaller nanoparticle size, the decay rates 
become similar above and below u s [see Fig. 3(b)]. This change in the frequency dependence 
is related to the stronger SP damping for u> > u s due to the d-hole scattering, as discussed 
above. Since this additional damping is reduced with decreasing temperature, the relaxation 
is faster above the resonance. This rather "nonlinear" relation between the time-evolution 
of the pump-probe signal and that of the temperature becomes even stronger for smaller 
sizes [see Fig. 3(c)]. In this case, the frequency dependence of the differential transmission 
decay below and above uj s is reversed. Note that a frequency dependence consistent with 
our calculations presented in Fig. 3(b) was, in fact, observed in the experiment of Ref. [T2| , 
shown in Fig. 3(d). 



VII. QUASIPARTICLE SCATTERING VIA SURFACE COLLECTIVE MODES 

Let us now turn to the electron scattering in the conduction band accompanied by the 
emission of surface collective modes. In the first order in the surface-induced potential, 
given by the second term in the rhs of Eq. (p2]) , the corresponding scattering rate can be 
obtained from the Matsubara self-energ; 
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1 4vre 2 \M LA f\ 2 

E - M w^-rSk ^' + <">• (30) 

where G c a = (iu — E^)~ l is the non-interacting Green function of the conduction electron. 
Here the matrix elements M%af are calculated with the one-electron wave functions V'aOO = 
Rni{ r )Yi m (r) . Since |a) and \a') are the initial and final states of the scattered electron, the 
main contribution to the Lth term of the angular momentum sum in Eq. ([3(]) will come from 
electron states with energy difference E a — E a i ~ u>l- Therefore, M^f can be expanded 
in terms of the small parameter Eq/\E^ — E c a ,\ ~ Eq/ul, where Eq = (2mR 2 )~ 1 is the 
characteristic confinement energy. The leading term can be obtained by using the following 
procedure. E We present M^f as 

M«< = (c,a^Y LM (r) M = (^J^j|l£X) ; (31) 

(E a - E a ,) 

where H = Hq + V(r) is the Hamiltonian of an electron in a nanoparticle with confining 
potential V(r) = V 9(r — R). Since [H, r L Y LM (r)] = — — V[r L Y LM (r)] ■ V, the numerator in 
Eq. (|3l|) contains a term proportional to the gradient of the confining potential, which peaks 
sharply at the surface. The corresponding contribution to the matrix element describes the 
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surface scattering of an electron making the L-pole transition between the states \c, a) 
and |c, a'), and gives the dominant term of the expansion. Thus, in the leading order in 
\E^ — E c a , we obtain 

m lm _ (c,aWr L Y LM (i)] ; W(r)|c,a') _ LR™ LM 



with tpimii m i = j c/rY i j! ra (r)Y^M(r)Yj/ m /(r). Note that, for L — 1, Eq. (j32|) becomes exact. For 
electron energies close to the Fermi level, ~ Ep, the radial quantum numbers are large, 
and the product VoR n i(R)R n 'i'( y R) can be evaluated by using semiclassical wave-functions. 
In the limit V — > oo, this product is given byi 2 J ' E c ^ l E^, u / R 3 , where E c nl = 7i 2 (n + l/2) 2 E 
is the electron eigenenergy for large n. Substituting this expression into Eq. fl3"2|) and then 
into Eq. (0), we obtain 

£«M = ~ E E E £4 t^t; 7#%ti S + *"). < 33 > 

P fa/ L n'V \ lL + ^ (^nZ - &n'l<) e L\^ ) 



with 



Civ = E l^k'l 2 = (2L+1 o K2r+1) /' dx^(x)P L (x)fl,(x), (34) 

where P(x) are Legendre polynomials; we used properties of the spherical harmonics in the 
derivation of Eq. (Q). For E° t ~ Pf, the typical angular momenta are large, I ~ k F R ^> 1, 
and one can use the large-/ asymptotics of Pj; for the low multipoles of interest, L -C I, 
the integral in Eq. fl34|) can be approximated by ^PT^'- After performing the Matsubara 
summation, we obtain for the imaginary part of the self-energy that determines the electron 
scattering rate 

lmKM = l iEgi{E)i EEl N { E-, ) + ! { E) 

where N(E) is the Bose distribution and gi(E) is the density of states of a conduction 
electron with angular momentum I, 




gi (E) = 2Y,6(E c nl -E)~-J—, (36) 



where we replaced the sum over n by an integral (the factor of 2 accounts for spin). 

Each term in the sum in the rhs of Eq. (|35D represents a channel of electron scattering 
mediated by a collective surface mode with angular momentum L. For low L, the difference 
between the energies of modes with successive values of L is larger than their widths, so 
that the different channels are well separated. Note that since all ujl are smaller than the 
frequency of the (undamped) bulk plasmon, one can replace €l(u) by in the integrand 

of Eq. fl55p for frequencies to ~ ojl- 

Consider now the L = 1 term in Eq. fl35|), which describes the SP-mediated scattering 
channel. The main contribution to the integral comes from the SP pole in ej" 1 (cu) = 3er 1 (c ( j), 
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where e s (u) is the same as in Eq. (|3]). To estimate the scattering rate, we approximate 
Ime^ 1 ^) by a Lorentzian, 

J a U%/u) 3 + (%(w) ^ u 2 u sl 



Ime c 1 \ 



10) 



[e'(uj) + 2e m } 2 + [wl/ufl + e'>(cu)} 2 e' d (co s ) + 2e m (^ - lj*)* + ^ 7 2 ' 

(37) 



where u s = lo\ = lo p / ^ e' d {uj s ) + 2e m and 7 = 7 S + LU s e d (u) s ) are the SP frequency and width, 
respectively. For typical widths 7 € w s , the integral in Eq. ( P^D can be easily evaluated, 
yielding 



24e 2 uj s E 2 E c a ^2m{u 



e d\ u s) + 2e m {uj — h/° — lo s ) 



Finally, using the relation e 2 kF[e' d {ui s ) + 2e m ] 1 = 3tiuj 2 /8E f , the SP-mediated scattering 
rate, 7|(££,) = — Im£^(£7£), takes the form 

Y e (E) = ( E ^ ± ) 1 ' 2 [1 " f(E ~ ".)]■ (39) 



LO s E F V 2?. 

Recalling that 2?o — {2mR 2 )~ 1 , we see that the scattering rate of a conduction electron 
is size-dependent: 7! oc 2?~ 4 . At 22 = 22p + o; s , the scattering rate jumps to the value 
97r(l + u s /Ef)Eq/u s , and then increases with energy as E z l 2 (for u s < £ F ). This should 
be contrasted with the usual (bulk) plasmon-mediated scattering, originating from the first 
term in Eq. (p2|), with the rate decreasing as E~ l l 2 above the onsetHI To estimate the size 
at which 7J becomes important, we should compare it with the Fermi liquid e-e scattering 
rateJEl 7 e (2£) = • F° r energies E ~ E F + uj S) the two rates become comparable 

for 

In the case of a Cu nanoparticle with uj s — 2.2 eV, we obtain k F R ~ 8, which corresponds to 
the radius R ~ 3 nm. At the same time, in this energy range, the width 7 | exceeds the mean 
level spacing 5, so that the energy spectrum is still continuous. The strong size dependence 
of 7! indicates that, although 7! increases with energy slower than 7 e , the SP-mediated 
scattering should dominate for nanometer-sized particles. Note that the size and energy 
dependences of scattering in different channels are similar. Therefore, the total scattering 
rate as a function of energy will represent a series of steps at the collective excitation energies 
E = ojl < uj p on top of a smooth energy increase. We expect that this effect could be 
observed experimentally in time-resolved two-photon photoemission measurements of size- 
selected cluster beams.E^I 



VIII. CONCLUSIONS 

To summarize, we have examined theoretically the role of size-dependent correlations in 
the electron relaxation in small metal particles. We identified a new mechanism of quasi- 
particle scattering, mediated by collective surface excitations, which originates from the 
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surface-induced dynamical screening of the e-e interactions. The behavior of the corre- 
sponding scattering rates with varying energy and temperature differs substantially from 
that in the bulk metal. In particular, in noble metal particles, the energy dependence of the 
o?-hole scattering rate was found similar to that of the absorption coefficient. This led us 
to a self-consistent scheme for the calculation of the absorption spectrum near the surface 
plasmon resonance. 

An important aspect of the SP-mediated scattering is its strong dependence on size. 
Our estimates show that it becomes comparable to the usual Fermi-liquid scattering in 
nanometer-sized particles. This size regime is, in fact, intermediate between "classical" 
particles with sizes larger than 10 nm, where the bulk-like behavior dominates, and very 
small clusters with only dozens of atoms, where the metallic properties are completely lost. 
Although the static properties of nanometer-sized particles are also size-dependent, the 
deviations from their bulk values do not change the qualitative features of the electron 
dynamics. In contrast, the size-dependent many-body effects, studied here, do affect the 
dynamics in a significant way during time scales comparable to the relaxation times. As 
we have shown, the SP-mediated interband scattering reveals itself in the transient pump- 
probe spectra. In particular, as the nanoparticle size decreases, the calculated time-resolved 
differential absorption develops a characteristic lineshape corresponding to a resonance blue- 
shift. At the same time, near the SP resonance, the scattering leads to a significant change 
in the frequency dependence of the relaxation time of the pump-probe signal, consistent 
with recent experiments. These results indicate the need for a systematic experimental 
studies of the size-dependence of the transient nonlinear optical response, as we approach 
the transition from boundary-constrained nanoparticles to molecular clusters. 

This work was supported by NSF CAREER award ECS-9703453, and, in part, by ONR 
Grant N00014-96-1-1042 and by Hitachi Ltd. 
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FIGURES 

FIG. 1. Calculated absorption spectra at positive time delays for nanoparticles with (a) R 
nm, (b) R = 2.5 nm, and (c) R = 1.2 nm. 



FIG. 2. Calculated differential transmission spectra at positive time delays for nanoparticles 
with (a) R = 5 nm, (b) R = 2.5 nm, and (c) R = 1.2 nm. 

FIG. 3. Temporal evolution of the differential transmission at frequencies close the SP reso- 
nance for nanoparticles with (a) R = 5 nm, (b) R = 2.5 nm, and (c) R = 1.2 nm. (d) Time-resolved 
pump-probe signal from Ref. GJI. 
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